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Abstract 

We discuss some special classes of canonical transformations of the extended phase 
space, which relate integrable systems with a common Lagrangian submanifold. Va- 
rious parametric forms of trajectories are associated with different integrals of motion, 
Lax equations, separated variables and action-angles variables. In this review we will 
discuss namely these induced transformations instead of the various parametric form 
of the geometric objects. 



1 Introduction 

Let us begin with a simple example. Consider an ellipse defined by the standard implicit 
equation 

— + = 1 
a? b 2 

One can represent this ellipse by the parametric equations 

x = asin(t), y = bcos(t), t £ [0,2tt]. (1.1) 

It is known, there are infinitely many parameterisations of a given curve. For instance, we 
can reparameterise an ellipse by using another parameter 

t= (a 2 + b 2 )t+ (a 2 - 6 2 ) sin(t). (1.2) 

Construction of different polynomial, rational and other parameterisations of the plane 
curves is a subject of classical algebraic geometry. 

In classical mechanics the same ellipse may be identified with integral trajectories of 
various integrable systems on the common phase space. In this case the parameter t is 
the time variable conjugated to some Hamilton function H. As an example, the first 
parametric form of the ellipse fll.lj ) is related to the two-dimensional oscillator, while the 



second parameterisation (1.2) may be associated to the Kepler model. 
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In this and several other unexpected situations in mathematics, dynamics is occasion- 
ally invading mathematical objects in which time is not present in the definition and 
yet the object can be endowed with various dynamical system structures, continuous or 
discrete. 

Thus, there is a problem of finding suitable parameterisation of curves, which are 
associated with various integrable systems. What do we know about this problem? 

We consider a mechanical system defined by some Lagrangian function C(q, q) or a 
Hamilton function H(p, q) on the 2n-dimensional phase space M with local coordinates 
{qj,Pj}^=i- According to Maupertuis' principle the extremals qj = "fj(t) of the action 
functional 

B 

C(q(t),q(t))dt, <7 =(<?!,... ,</n) (1.3) 

A 

coincide with the extremals of the reduced action functional Sq on the fixed energy surface 

Q 2n - 1 =(H(p,q) = E). (1.4) 

Recall that the Maupertuis' principle was at first enunciated in 1744 [2(1, IS, [10]. The 
modern interpretation of the Maupertuis' principle may be found in [p|, P2,p||. The reduced 
action Sq = J pdq is independent of any evolution parameter. Moreover, even its initial 
parameter t and the corresponding Hamilton function cannot be restored from the reduced 
action problem Q. Nevertheless, solutions of the corresponding variational problem are 
the initial trajectories qj = Jj(t) in the common nonparametric form §. For instance, 
trajectories of the Kepler system are conic sections ar^ 1 = 1 + bcoscj), which may be 
ellipses, hyperbolas or parabolas at E < 0, E > and E = 0, respectively. 

So, any integral trajectory may be parameterised by using another parameter t, such 
that 

ij = 7i(t) = te[A,B], te [C, D]. 

According to Maupertuis' principle all the initial trajectories have one nonparametric form 
on the surface Q 2n_1 . Therefore, for all these trajectories we could introduce common 
parametric form qj = f3j(t) and a new local Hamilton function H defined on Q 2n_1 . Two 
main problems are how to find the new parametric form of the initial trajectories and how 
to obtain the new Hamilton function defined on the whole phase space A4. 
Consider the corresponding Hamilton-Jacobi equation 

dS fdS_ dS \ 
dt \dqi"' ' dq n ) 

In the invariant geometric Hamilton-Jacobi theory [^4| any hyperplane Q 2n_1 in Ai may 
be called the Hamilton-Jacobi equation. Its solution is an ra-dimensional Lagrangian sub- 
manifold C( n ) in A4, such that C Q 2n_1 . By definition, a Lagrangian submanifold is 
one where the symplectic form Q vanishes when restricted to it, i.e. Q\ c = 0. This defi- 
nition is completely invariant with respect to change of local coordinates and parametric 
representations of the Lagrangian submanifold As above, we could consider various 

parameterisations of a given Lagrangian submanifold. Each new parametric form yields a 
new Hamiltonian system related to the same geometric surface. 
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Below we consider integrable systems on M with n integrals in involution. According 
to the Liouville theorem Q for any integrable system the corresponding n-dimensional 
Lagrangian submanifold depends at least of n- arbitrary constants. So integrability is a 
geometric property and it does not depend on the choice of the parameterisation of the 
Lagrangian surface. Starting with a known Lagrangian surface of some integrable system 
we can try to get new integrable models by using various parametric forms of this surface. 
In this case, we can expect that the initial and resulting integrable systems have a lot of 
common properties. The main problem is how to find different parameterisations and the 
corresponding sets of integrals of motion defined on the whole phase space Ai . Generally 
there is no rule for how to proceed. Each case is different. 

Usually the Lagrangian submanifold depends on m > n arbitrary constants. The n 
constants a±, . . . ,a n are identified with the values of integrals of motion Ij = otj [44], 
while the remaining m — n constants a±, . . . , a m - n are free parameters. In this case we 
have a freedom related to the choice of n integrals of motion from the n+m initial constants 
of motion. This freedom permits us to associate a family of integrable systems with one 
Lagrangian surface. In this review we discuss a special class of different parametric forms 
of a given surface, which is associated with mutual permutations of energy E and arbitrary 
parameter a&. The corresponding reparameterisation of we will call the generalized 
Kepler change of time. 

The aim of this paper is to bring together some old and new examples of integrable 
systems related with the various parametric forms of the one Lagrangian surface. The 
passage from a given parameterisation to the another one gives rise to the transformations 
of all the properties of integrable systems, such as integrals of motion, Lax equations, 
separated variables and the action-angles variables. In this paper we discuss namely these 
induced transformations instead of the various parametric form of the geometric objects. 

The initial symplectic form Q is equal to zero on the Lagrangian submanifold . We 
can suppose, that the space M may be equipped with another form f2, which is equal 
to zero on the same surface C^ n \ In this case for a given Lagrangian surface time t and 
symplectic structure may be changed simultaneously. Moreover, we can try to embed 
this n-dimensional surface into the various phase spaces. For instance, the ra-dimensional 
Kepler problem may be identified with the geodesic flow on the sphere S n [17, 21]. The 
Kolosoff transformation maps the Kowalevski top into the two-dimensional Stackel sys- 
tem [f[l[] . The same top may be related with the geodesic motion on SO (4) ||] or with 
the Neumann system on the sphere S 2 [11]. Here we will not discuss such composed 
transformations of time and phase space. 

By embedding a Lagrangian submanifold into the infinite-dimensional phase space, 
we can identify with an invariant manifold of some hierarchy of nonlinear evolution 
equations. Such finite dimensional manifolds are invariant with respect to the action of 
all flows of the hierarchy and they are naturally expected to be integrable since all the 
flows of hierarchy commute on them. For instance, given Lagrangian submanifolds may 
be realised via stationary or restricted flows. Here we will not discuss relation of finite- 
dimensional integrable systems with soliton equations and restrict ourselves to finding new 
parameterisations of the known trajectories only. 

Below we will consider many well known mechanical systems, such as the Stackel sys- 
tems, Toda lattices, Henon-Heiles and Holt systems, integrable systems with the quartic 
potential and the Goryachev-Chaplygin top. The results we present are, for the most 
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part, not new and we do not provide detailed proofs (these can be found in the papers 
cited). What may be new and interesting is an exposition of canonical transformations 
of the extended phase space as different parametric forms of integrable geometric objects 
and the action of these transformations on the properties of integrable systems. 



2 The Maupertuis— Jacobi transformations 

Let M be an n-dimensional Riemannian manifold with the metric gij. On the cotangent 
bundle A4 = T*M consider a Hamiltonian system with the natural Hamilton function 
H(p,q) 

n 

H(p,q) = T(p,q) + V{q) = ^gMjHPj + V(q). (2.1) 

On the smooth submanifold Q 2n_1 (|1.4|), integral trajectories of the Hamiltonian vec- 
tor field £ = sgrad H(p, q) coincide with integral trajectories of another vector field £ = 
sgrad H{p, q), where the new Hamilton function is given by 

H{p, q) = Tip, q)=J2 E l iv\ q) P*Pr ( 2 - 2 ) 

U 

Integral trajectories have two different parametric forms qj = 7j(i) = (3j(t) on the surface 
q2?i-i on iy_ However, the resulting Hamilton function describes geodesic motion and, 
therefore, on the whole phase space we can determine the so-called Maupertuis transfor- 
mation 

which relates the initial hamiltonian vector field £ on A4 with the other hamiltonian vector 
field £ defined on the same phase space M. ||. 

If t be the time along trajectories of the new vector field £, then the Maupertuis mapping 



gives rise to so-called Jacobi transformations [10, 16, 27] of the Hamilton function ( p.2[ ) 
and of the time variable 

dt = (E - V{q))dt. (2.3) 

This Jacobi transformation explicitly describes new parameterisation of the common in- 
tegral trajectories q(t) = q(t) and determines the new Hamilton function H(p,q) (|2.2| ). 

The Maupertuis transformation maps any integrable system with a natural Hamilton 
function H (p, q) into the other integrable system on the same phase space Ai. Namely this 
property has been used for the search of the new integrable systems (see references within 
[19, 1C, l|, |28|). The Maupertuis principle for integrable systems with a nonnatural 



Hamilton function is discussed in 

The property of integrability is independent on the choice of parametric form of tra- 
jectories. However, some criteria of integrability drastically depend upon the parameteri- 
sation. For instance, the method of singularity analysis associates integrability with the 
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Kowalevski-Painleve property, i.e. the only singularities of the solutions of the equations 
of motion are movable poles (t — to) _m i n the complex t-plane |l], |23]]. There exist cases of 
integrable systems with the rational integrals of motion, whose analytic structure permits 
solutions with algebraic singularities of the type (t — to)~ m ^ k , (k being a positive integer 
larger than one). These systems satisfy to the so-called "weak" Painleve property (see 
review ||). 

We have some examples of integrable systems related with one Lagrangian submani- 
fold [23, 39], which satisfy to the usual Kowalevski-Painleve property and the "weak" 
Painleve property, respectively. These systems are related with the different parametrisa- 
tions of one geometric object. So we can see that a change of parametric form leads to a 
change of the Kowalevski-Painleve criteria of integrability. 



3 Canonical transformation of the extended phase space 

To find new parameterisation of the known integral trajectories or of the Lagrangian sur- 
faces one has to introduce a new parameter t and the corresponding Hamilton function H. 
Thus, to describe explicitly the following mapping 

t^t, H{p,q)^ H(p,q) (3.1) 

we extend initial phase space M by adding to it the new coordinate q n +i = t with the 
corresponding momentum = —H. The resulting 2n + 2-dimensional space Me 

PI |34j is the so-called extended phase space of the hamiltonian system. We emphasize 
that H{p,q) is the Hamilton function on M, but H is an independent variable in the 
space Me- The energy E is a fixed value of the variable H or the function H(p, q). 

To describe evolution on the extended phase space we introduce the generalized 



Hamilton function [16, 34] 



H(p!,.. . ,Pn+i]qi ■ ■ ■ ,Qn+l) = H(p,q) - H. (3.2) 
The Hamilton equations for the variables q n +i = t and p n +i = —H are 

dt _ 1 dH _ o 
dr dr 

Here r is a generalized time (parameter) associated to the generalized Hamilton func- 
tion Ji. The time variable t is a cyclic coordinate and the conjugated momentum is a 
constant of motion. Other 2n equations coincide with the initial Hamilton equations on 
the zero-valued energy surface 

H(p,q)=H(p,q)-H = 0. (3.3) 

Thus, our initial hamiltonian system on A4 may be immersed into the hamiltonian system 
on Me ■ Using this immersion and canonical transformations of the extended phase space 



Me [16, e3 we introduce transformations 

(H, t) canonical transformations (H, t) 
M E — > M E 
T I 
M M 

H(p, q) H(p, q) 
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which map an initial Hamilton function H (p, q) \— > H(p, q) into another Hamilton function 
defined on the same phase space. Of course the similar mapping t t— > t permits us to 
describe new parameterisation of the corresponding integral trajectories. 

Note two different classical definitions of the canonical transformations are known Q. 

1. Canonical transformations preserve the canonical form of the Hamilton- Jacobi equa- 
tions. 

n 

2. Canonical transformations preserve the differential 2-form, £1=^2 dpif\dqi, on M. 

i=l 

For instance the first definition is used in textbooks on variational principles of classical 
mechanics [1C, [34], [|. The second definition was later introduced to consider geometry 
of the phase space 44 1 . 

Below we use the first definition of canonical transformations because the Maupertuis- 
Jacobi transformation and the Kepler change of the time preserve the canonical form of 
the Hamilton-Jacobi equation, but it retains the corresponding differential 2-form $7 = 

n+l 

£ dpi A dqi on the level Q 2n_1 (0) only §]. 
i=i 

We introduce general canonical transformations of the extended phase space M. e 



(3.4) 



t i— ► t, dt = v(p, q)dt, 
H^H, H = v{p,q)- l H. 
which change the initial equations of motion 

dqi -i ( s f dqi ^dv\ dpi _ x ( dpi ~dv 

but preserve the canonical form of the Hamilton-Jacobi equation 

— + H = 0, where S= (pdq-Hdt). (3.5) 

and retain the corresponding zero-energy surfaces ( |3.3[ ) at v(p,q) ^ 

H(p,q) =v(p,q)- 1 H(p,q) =0. 

Zeroes of the function v(p, q) determine the behavior of the system with respect to the 
inversion of time Q. Here we do not consider this problem in detail. 

The Maupertuis-Jacobi mapping (|2.2| ), (|2.3| ) may be rewritten as such a canonical 
transformation ( [SOI ) of the extended phase space 

T{p, q) i ^ Tip, q) = v(p, g) _1 T(p, q), v{p, q) = E - V(q), 

which maps an initial geodesic flow into another geodesic flow. This map preserves inte- 
grability, if the function v(p, q) is constructed by any potential V(q), which may be added 
to the initial kinetic energy T(p,q) without loss of integrability [10]. 

In contrast with the Maupertuis-Jacobi transformations, even if the general canonical 
transformation of A4e fl3.4| ) retains integrability, we have no general method to construct 



new integrals of motion starting with initial ones. To solve this problem in 39, 41, 42, 
43], we used some analogies with the one known example of such transformations due to 
Kepler. 
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4 The Kepler change of the time 

We begin with brief description of the Kepler change of the time [17]. We commence with 
two-dimensional oscillator defined by the Hamilton function 

H osc (p, q)=vl+pl + a (qj + qf) + b, a, b € R. 

For this system the Kepler canonical transformation (|3.4| ) of Me with the function 

v(p,q)=qf+ql (4.1) 

preserves integrability. After change of the time (^J) and the point canonical transforma- 
tion to other variables 

x = qm, y = {q\- ql) /2 (4.2) 

integral trajectories of the oscillator come to be trajectories of the Kepler problem defined 
by 

H kep i{p, x) = =P X +P y + —=== + a. (4.3) 

<?i + ?2 2y ar + y z 

Various parametric and nonparametric forms of the common trajectories are discussed 
in [17, ^, |j. Coincidence of the integral trajectories may be regarded as a local result, 
whereas the corresponding canonical transformation of the extended phase space preserves 
integrability in the whole initial phase space. 

As for the Maupertuis- J acobi transformation, the function v (p, q) ( |4.1| ) in the Kepler 



transformation (L3) could be identified with the oscillator potential V(q). Below we prove 
that it is a simple coincidence. Nevertheless canonical transformations of the type (3.4) 
have been called the coupling constant metamorphoses in [jl^, ^] because of this casual 
coincidence. 

The Kepler change of time has been generalized by Liouville [18]. The Liouville inte- 
grable systems are systems with the natural Hamilton function 

H(p,q)=f + V, 

where the kinetic and potential energies are given by 



'v, 

i=l i=l i=l 



T = v-\q)Y,aiPl V = v-\q)J2u i , v(q) = £\ 



The functions a^, Vi and Ui depends only on the variable qi. 
For the Liouville system the following quantities 

Ij = ajPj + Uj - Hvj, j = 1, . . . , n, 

are integrals of motion in involution and ^ Ij = 0. Thus, we have n quadratic integrals 
of motion including the Hamilton function and consequently the Liouville systems are 
completely integrable. Note the two-dimensional oscillator and the Kepler model belong 
to the Liouville family of integrable systems. 
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We recall how equations of motion were integrated in quadratures by Liouville [18] 
From Ij = otj, one obtains a system of differential equations 

dqj dt 

7 = 1, ... ,n. 



a 3 [ a j + Ev j- U j 



v{q) 



Here H = E and time t is associated with the Hamilton function H. Choosing a new time 
variable t according to 

dt = v~ l (q)t 

we come down to the system of equations 



13 



dt. 



aj [aj + Evj - Uj 



It allows one to find integral trajectories qj = Jj(t). After that the new parameter t may 
be expressed in terms of the initial time t by the quadrature 

t = J v(qi(r),... ,q n (r))dT. (4.4) 

This transformation is related to a new parametric form of the same trajectories qj = /3j(t). 

In fact, Liouville tacitly used canonical transformation of the extended phase space ( |3.4[) 
and considered a new integrable Hamilton function H = v{q) H instead of the initial one. 
After integration of equations of motion for the new system we have change parametric 
form of the trajectories in order to describe solution of the initial problem. For instance 
in parametric form of the ellipse the Kepler time ( |1.2| ) is explicitly the Liouville quadra- 
ture flop. 



Note canonical transformations of the extended phase space ( |3.4|) have a natural coun- 
terpart in quantum mechanics. Namely, it is known that the standard eigenvalue problem 
of the Hamiltonian operator H(p,q), 

H(p, q)V = (H + aV + b)^> = m, 

may be associated with the eigenvalue problem of the charge operator a 

H(p, q)V = (h + (b- E)V~ l \ $ = -a*. 



In quantum mechanics such a duality of the two eigenvalue problems has been used 
by Schrodinger and many other [Pq] . Canonical transformation of the extended phase 
space ( |3.4D is an analogue of this duality. It is interesting that for the first time this duali- 
ty has been studied for the quantum Kepler problem. In the Birman-Schwinger formalism 
function v(q) is called a "sandwich" potential p5fl . Recall that in the Birman-Schwinger 
formalism we can estimate the spectrum and eigenfunctions of the one Hamiltonian H, by 
using the known spectrum and eigenfunctions of the dual Hamiltonian H. Moreover, for 
some quantum models it is a single known way to find solutions of the initial Schrodinger 
equation. Below we briefly discuss a similar property for the quantum Toda lattice. 
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5 The generalized Kepler change of time 

The Maupertuis-Jacobi mapping is traditionally used for the search of new integrable 
systems. The Kepler change of time and the Liouville reparameterisation of trajectories 
have been used for integration of equations of motion. In this Section we propose some 
generalisations of the Kepler-Liouville results, which may be useful for the search of new 
integrable systems as well. 

All the Liouville systems are particular case of Stackel integrable systems. Therefore 
let us briefly recall some necessary facts about Stackel systems ^M- The nondegenerate 
n x n Stackel matrix S, its j-th column of which Skj, depends only on qj 

detS/0, ^ = 0, j^m 
oq m 

defines the set of functionally independent integrals of motion, {Jfe}^ =1 , where 

h = Yl c i k (P 2 J + U J fe")) » c jk = ^~ g , (5-1) 

3=1 

which are quadratic in the momenta. Here C = [cjk] denotes the inverse matrix of S 
and S fej is the cofactor of the element Sjy. 



Proposition 1 [38]. If the two Stackel matrices S and S be distinguished by the m-th 
row only, i.e. 

Skj = Sfcj, k ^ m, 

the corresponding Hamilton functions I m and I m ( \5.1[ ) with a common set of potentials Uj 
are related by canonical transformations of the extended phase space M e 

I m i-> Im = v~ 1 (q)Im(p, <?), dt m = v(q)dt m , (5.2) 

where 

detS(gi, ... ,q n ) (K ^ 

v(q u ... ,q n ) = — — -. (5.3) 

detS(gi, ... ,q n ) 

The proposed generalization of the Kepler transformation maps one integrable Stackel 
system into another integrable Stackel system. For instance, the Stackel matrices for the 
oscillator and the Kepler problem are 



It is obvious that the Kepler change of time ( [4.3[ ) coincides with the proposed map- 
ping (|5~^). 

We return to the Kepler transformation ( |4.3j ) of Me- After permutation of coordinates 
and momenta (q\ t 2 ^ Pi, 2) the Hamilton function for the Kepler problem 



Pi + V\ + a^ 1 {q\ + qfj + a~ l b _ H, 



Hkepl — a 2 1 2 — a 



osc 



p{ + pi Hfrcc 
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becomes a ratio of the Hamilton function H osc for the oscillator and the Hamilton function 
Hf TCC for the free motion. 

So for any two integrable systems the ratio of their Hamilton functions could be the 
Hamilton function of a third integrable system on the same phase space. The main re- 
maining problem is a search of a complete set of integrals of motion. 

We consider two integrable hamiltonian systems on the common phase space Ai . These 
systems are defined by the two sets of independent integrals of motion {Ij}™ =1 , and { Jj}™ =l 
in involution, i.e. 



{Ij,Ik} = and { Jj, J k } = 0, j, k = 1, 



n. 



Introduce the antisymmetric matrix /C = (I ® J), which is the inner product of the two 
independent vectors of integrals I and J in M n . Any column or row of this matrix defines 
a set of n — 1 independent functions 

K-%j = (I ® J)jj = li Jj ~ Ij J%i 2, j = 1, . . . , n. 



Proposition 2 |41|| . If all the differences of integrals of motion (Ij — Jj) with the common 
index j = 1, . . . ,n are in involution, i.e. 

{Ij - Jj,I k - J k } = , j, k = 1, . . . ,n, (5.5) 

then the ratio of integrals 

K m = 1 -f- (5.6) 
and n — 1 functions Kj , j ^ m 

Kj = — — = '—j— 1 = K m Jj - Ij, m / j = 1, . . . , n (5.7) 

are integrals of motion for new integrable system on the same phase space. 



Thus the mapping ( |5.6D defines a canonical transformation (3.4) of the extended phase 
space, which preserves the property of integrability. To apply this transformation we have 
to find two integrable systems satisfying condition (p^). 

We consider a pair of the Stackel systems with a common Stackel matrix S and with 
different potentials Uj. Namely, in addition to the system with integrals {I k } ( |5.1[ ), we 
introduce the second integrable system with the similar integrals of motion, 

n 

J k = c ik {P 2 j + Wj( qj )) , k = 1, . . . , n. (5.8) 
i=i 

At least one potential Uj(qj) has to be functionally independent of the corresponding 
potential Wj{qj). 



Proposition 3 |41 |. Any two integrable systems defined by the same Stackel matrix S 
and by functionally independent potentials Uj{qj) and Wj(qj) satisfy the necessary con- 
dition /f5. of the previous proposition. Thus, the ratio of the two Stackel integrable 



Hamiltonians defines new integrable system 

(I m ,Jm) K m = v(p,q)~ 1 I m = (5.9) 
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It is obvious that all the integrals Ik and Jk differe by the potential part 

2 

(4 - J k ) = ^2 c ik PMi) - w j(u)} 

3=1 

depending on the coordinates q only. Thus systems with a common Stackel matrix S 
satisfy condition (5.5). 

The Hamilton function ( |5.9D has the following form 



H(p,q)=K v 



E c jm p] + Uj(qj) - P n 



E9m Pj + Wjiqj) 
i=i L 



Pm € 



(5.10) 



This Hamiltonian H (p, q) is a rational function in the momenta, but next one can try to 
use canonical transformations to simplify it. Occasionally, one obtains again a natural 



type of Hamilton function. For instance, according to [41], integrable systems with the 
following Hamilton functions 



Hi =p k x p k y + a{xy) * 



k 

+ 1 



a,k G 



(5.11) 



H n =p k +p k + a(xy) fc+i, 



belong to the proposed family of generalized Stackel systems. At k = 1 the first Hamilto- 
nian coincides with the Hamiltonian of the Kepler problem. At k = 2 the second integrable 
Hamiltonian has been found by Fokas and Lagerstrom (see references within |23|]). In this 
case both initial systems satisfy to Kowalevski-Painleve criterion, whereas the resulting 



Fokas-Lagerstrom system admits asymptotic solutions with fractional powers in t [23|. 



6 Properties of the change of time for the Stackel systems 

For the Stackel family of integrable systems we proposed two different examples, ( |5.2| 



and (|5.9|), of the canonical transformations (3.4) of the extended phase space Me- Now 
we consider some properties of these transformations. 

Recall that the common level surface of the Stackel integrals Ij 

{zeM 211 :I l (z)=a i , i = 1, . . . ,n} (6.1) 



is diffeomorphic to the n-dimensional real torus. We can immediately construct the one- 
dimensional separated equations 

2 n 



Pi 



( dS \ 

\~dqj) = = ^ aiSi ^ ~ U i( q i> a )' ak G 



(6.2) 



Here Sq is a reduced action functional [33]. Integral trajectories qj(t, a\, . . . ,a n ) are 
determined from the following equations 



E 

j'=i' 



Skj{qj)dqj 



Pk, Pi = t. 



(6.3) 
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In fact the polynomial P(X) and the contour of integration depend upon the values aj of 
the integrals of motion, which are dropped in the notation. 

For the rational entries of S and rational potentials Uj(qj) the Riemann surfaces ( |G.2| ) 
are isomorphic to the hyperelliptic curves 

2g+l 

Cj : = P(Xj) = ]T a k \ k 3 . (6.4) 

k=l 

Considered together these curves determine the n-dimensional Lagrangian submanifold in 
the phase space M = M 2n 

: Ci(pi, qi) x C 2 {p 2 , 92 ) x • • • x C n (p n , q n ), (6.5) 

which is decomposed into plane curves. Applying Arnold's method [^| we find that the 
action variables have the form 

5j = I P(Xj) d\j. (6.6) 
J Aj 

The Abel transformation linearizes the equations of motion on the Lagrangian subman- 
ifold (j( n ) in terms of abelian differentials of the first kind on the corresponding spectral 
curves [37, 38]. 



We consider a pair of the Stackel systems related by the first generalization of the 
Kepler mapping (5.2) such that the potentials Uj(X) = ^a k X k depend on arbitrary pa- 



rameters a k . According to initial and resulting integrable systems are associated with 
algebraic hyperelliptic curves (|6.4j) C and C described by 



C: n 2 =J2 a J Xi + a mX m + EX k + a i X ^ 

_ _ (6-7) 

C: = ^ j a j X : > + EX m + a k X k + ^aiX\ 

Here the n coefficients {aj} and {5j} are values of the integrals of motion such that E = a± 
and E = ol\. The other coefficients aj G M are arbitrary parameters (charges), which 
define the potential part of the Hamilton function H(p,q) = T(p) + V(q,a). Canonical 
transformations of Me give rise to mutual permutation of the energy E and one of the 
parameters a k (charge). 

The initial and resulting Riemann surfaces are topologically equivalent. One can prove 
that the corresponding integrable systems are topologically equivalent too. Moreover, 
initial curves coincide with the resulting curves at the special values of integrals of motion 

E = a m , E = a k , aj = aj, j = 2,... ,n. 

In this case we have two different parametric forms of the common Lagrangian submani- 
fold CW (H) depending on n values of integrals aj and constants a k - On the correspond- 
ing submanifolds A4 a and M. a (13) integral trajectories of the initial system qj(t) fl6.3|) , 
coincide with integral trajectories of the resulting system qj(t). In the neighbourhood of 
the intersection of these submanifolds we can introduce the common set of the action- 
angle variables ( |6.6[) . In this region the Kepler transformation (^^) retains the differential 
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2-form Q in Me and the function v(p, q) = v(s) depends on the common action variables 
only. 

We consider a pair of the Stackel systems related by the second generalization of the 



Kepler mapping (5.9). According to |^l[ initial and resulting integrable systems are asso- 
ciated with algebraic hyperelliptic curves (]6.4|) C and C of the form 



C: fi 2 = Y^a i X i + Y, a i Xi ^ C: ^ = ^b k X h +^a m X m . 

Here the resulting coefficients are rational functions of the initial ones [41]. So canonical 
transformations of Me ( |5.9| ) give rise to transformations of the modulus of the curves 
only. 

As mentioned above the corresponding Riemann surfaces and integrable systems are 
topologically equivalent. At the special choice of integrals {a, a} and parameters of po- 
tentials {a,, b} integral trajectories of the initial system coincide with trajectories of the 
resulting system. In a small region of M this generalization of the Kepler transforma- 
tion ( |5.9D preserves the differential 2-form O in Me and the function v(p,q) depends on 
the action variables only. 

In an attempt to understand the origin of the conservation of integrability by canonical 
transformation of Me ( P-4| ) one can attempt to rewrite equations of motion in the Lax 
form 

{H(p,q),L(X)} = [L(X),A(X)}. 



For the some classes of Stackel systems the Lax matrices have been constructed in j37|. |38|, 
40]. 

We consider some examples only. One of the simplest Lax matrix L(X) for the one- 
dimensional Stackel systems is given by 



P 



X 



L(X) 



X-q 



(6i 



AIN 



-P 



Here <fi(X) = ^ (f>k^ k is a parametric function of the spectral parameter A and the elements, 
Nmat, are the linear combinations of the Taylor projections 



[ z \n 



£ 

k=— oc 



ZkX 



N 



N 

^2*kX\ 

k=0 



(6.9) 



or the Laurent projections [pq |. The coefficients <pk of the function ^(A) are the constants 
of motion, which may be parameters or integrals of motion. So at arbitrary M, N the 
family of the Lax matrices L l p(X) may be associated with a web of algebraic curves instead 
of one concrete curve. 

For instance, we present some one-dimensional Hamilton functions, the functions 4>(X) 
and the corresponding Lax matrices associated with the first Taylor projection [z]oi- An 
application of the pure numeric function <f>(X) yields 

k2 



H 



L 



p 2 + aq 2 + bq, 



p 

-a(X + q) 



X 



-P 



-aX z - bX, 



A 







(6.10) 
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The corresponding spectral curve is given by 



C: n 2 = a\ 2 + b\- H. 

Note that in the spectral curve we can substitute integrals of motion, while in the corre- 
sponding Lagrangian submanifold we have to substitute the values of integrals of motion. 
From ( |6,7D the first possible change of this curve looks like 

C i-> C : /u 2 = aA 2 + (6 - H)X + c. 
The associated canonical transformation of the extended phase space, 

H = v~ 1 (H + c) =q' 1 (H + c), 
changes the Lax matrices by the following rule 





-aX 2 - {b- H)X, L = L + H 



1 



,4 



-\q)A. 



(6.11) 



The second possible change of the curve 

C i — * C : /u 2 = (a - H)X 2 + bX + c 
may be related to the other canonical transformation of the extended phase space, 

H = v~ 1 (H + c) =q~ 2 (H + c), 
such that 

bX, 



4> = -(a - H)X 2 
L = L + H 





l + q 



A 



-Hi) 



A + H 





1 



(6.12) 



It is not hard to check that the Poisson bracket relations for all these Lax matrices L(X), 
L(X) and L(X) are closed into the linear r-matrix algebra. At the first case the r-matrix 
is a constant matrix r = 11/ (A — /i), whereas the second and third r-matrices depend on 
dynamical variables. Here II is a permutation of auxiliary spaces Sj. 

We turn now to the original change of time in Kepler problem (JD^) . The Lax matrices 
for the two-dimensional oscillator, 



C(X) 






^2(A,p 2 ,g2) 



A(X) 




(6.13) 



may be constructed from two independent 2x2 blocks Lj(X) and Aj(X) of ( 6.1CQ . The 
corresponding spectral curve C = C\ x C2 is a product of two hyperelliptic curves. 

The Kepler mapping ( |4.3] ) gives rise to the following transformation of the Lax matrices 

/ \ 

X + qi 



V H92 0/ 



£(A) = £(A) - H kcpl 



A{X)=v{p,q)- 



A{X) - H] 



kepi 



/ \ 

10 



V 1 / 



(6.14) 
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The spectral curve C remains a product of two new hyperelliptic curves. The initial 
oscillator and the resulting Kepler model are separable in the common variables, which 
lie on these curves. Note that these separated variables are cartesian coordinates for the 
oscillator and parabolic coordinates for the Kepler problem. 

Similar transformations of the Lax matrices ( |6.11 ), ( |6.12[ ), ( |6.14 ) may be proposed for 
the other two-dimensional Stackel systems separable in cartesian or parabolic coordinates. 
For the two-dimensional Stackel systems separable in elliptic or polar coordinates trans- 
formation of the Lax matrices has a more complicated form |38[ |. These transformation 
may be constructed by using two different outer automorphisms of infinite-dimensional 
representations of underlying si (2) algebra proposed in [ 40 1 . 

In the next Sections we consider similar transformations of the Lax matrices ( |6.11| ) and 
the spectral curves ( |6.7D for integrable systems associated with non-hyperelliptic algebraic 
curves and for the non-Stackel integrable systems. 



7 On integrable systems with quartic potential 



According to [23, |3Sj, f43[, canonical transformations of the extended phase space relate 
three integrable cases of Henon-Heiles systems with the three integrable cases of the Holt 
systems. One of these systems admits a 2 x 2 Lax matrix and the corresponding spectral 
curve is an hyperelliptic algebraic curve. Another two systems possess 3x3 Lax matrices 
and the corresponding spectral curves are the trigonal algebraic curves = P(X). As for 
the Stackel systems transformations of these 3x3 Lax matrices are shifts of their entries 
by the element of the extended phase space ( |6.11 ). 

We consider two integrable systems with quartic potential for which 4x4 Lax matrices 
were constructed in B by applying relations with stationary flows of some known inte- 
grable PDEs. The first system belongs to the Stackel family and is separable in cartesian 
coordinates. Its Hamilton function is 



H(p, q)=p 2 l+ p 2 2 + ^ (qf + §q\q\ + g£) 



(7.1) 



As above we could construct some block 4x4 Lax matrices (6.13) for this system in carte- 
sian separated variables. The corresponding spectral curve is a product of hyperelliptic 
curves. The two pairs of the separated variables lie on these two curves. 

Another 4x4 Lax matrix for the same system may be obtained from the Lax repre- 
sentation for the Hirota-Satsuma coupled KdV system Gj. This Lax matrix 



L(A) 



/ 


-PlQl 


8? 





1 \ 






Pi Qi 


2 










A 


-P2Q2 


4 


V 


A 2 





\-pI 


P2Q2 ) 



does not have a pure block-diagonal structure and the corresponding spectral curve /U 4 = 
P(A) is not a product of two hyperelliptic curves. The second Lax matrix A(X) may be 
found in Q. 
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+ vtr 1 (h - b) 



For this Stackel system canonical transformation of the extended phase space (5.2) 
H 

gives rise to the shift of the first Lax matrix and the rescaling of the second Lax matrix 

/ \ 

10 



\ A / 



L = L + H 



A = v~ x A. 



In contrast with the Stackel systems transformation of the corresponding spectral 

curves acts on the both side of the equation defining the curve 

C: fi 4 = A 3 - HX 2 + JA, 

(7.2) 



C 



iu 4 - 2fi 2 XH = A 3 



b) A + JA. 



Note that after canonical transformations of the other variables (p, q) the new Hamilton 
function H may be rewritten in the natural form 

~ 2 2 , x 2 + 2y 2 -b 

H = Px + Pv H / ■ 

y 2^x 2 + y 2 

The second integrable system with the quartic potential is the non-Stackel system 
defined by the following Hamilton function 

H = v\ + V\ ~ I {qf + Qqh 2 2 + 8qt) ■ (7.3) 

This system is separable after a so-called quasi-point canonical transformation p3]. 
The second system possesses 4x4 Lax matrix 

/ 



L 



^ + qiPl 



-q\ 



2q 2 



q ^+p\ + qiq2Pl + h ~ q -t-qm 



-2q 2 X 
{ A (-p2 +q 2 2 + 4 



2X 




P2 + q 2 2 + 9 4 
- q -^- + qiPl 



2 




-q\ 



qI f + Pi - qmpi + I ^f-qiPi ) 



which was obtained using a gauge transformation of the Hirota-Satsuma coupled KdV 
system §. 

For this system we propose the canonical transformation of the extended phase space 
given by 

H={q 2 +Aq 2 )- 1 (H-b), 

which preserves integrability. As for the Stackel systems, transformation of the Lax ma- 
trices retains a shift of the first Lax matrix and rescaling of the second Lax matrix 



L = L + 2H 



( \ 

10 



V A / 



A = v~ 1 A. 
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The spectral curves are changed by the rule 
C: fi 4 = A 3 + 4HX 2 + JA, 

C: /J, 4 - 8fi 2 \H = A 3 + 4 ( b - AH 2 ] X 2 + JA. ^'^ 



Recall that for Stackel systems with quadratic integrals of motion canonical transforma- 
tions of the extended phase space (|5.2| ), ( fTsI ) preserve separated variables lying on the 
common hyperelliptic curve (|6.7|). 

For the Henon-Heiles systems and systems with quartic potentials transformations of 
the trigonal spectral curve /x 3 = -P(A) |?9| and the curve // = P(X) (7^2), (|7.4| ) have a more 
complicated form. In this case separated variables for the initial and resulting systems are 
different [j39||. It means that the proposed canonical transformation of the extended phase 



space preserve integrability, but seriously changes other properties of the systems. 



8 The Toda lattices 

Before proceeding further, it is useful to recall some known facts about the Toda lattices 
(all details may be founded in the review p6[ ). 

Let q be a real, split, simple Lie algebra of rankg = n and let K(-, •) be its Killing form 
and P be a system of simple roots, we identify the phase space with the coadjoint algebra 
A4 ~ q* r = q+ Here q + is a Borel subalgebra and g_ is the opposite nilpotent 
subalgebra of q. 

The Lax matrices for the Toda lattices are 



£ = ^2 pihi + ■ exp K ( a > ^ • e -a + 

a a e ai 

i=l aSP a£P /g 

A = — y~] exp K(a, q) ■ e_ a , a a G R. 

The Hamilton function is given by 

H{p : q) = \k(C,C) = l -K{p,p) + Y, ^ a(q) - (8-2) 

Recall that in a shifted version of the Adler-Kostant-Symes scheme in order to construct 
the Toda orbit ( |8.1|) in A4 we have to translate a dynamical orbit living in q*^_ by adding 
to it a constant vector e = ^ a a e a from the remaining part of A4. This vector e has to 
be a character and has to be a constant. 

We replace now the phase space M. by the extended phase space M. e ■ Roughly speak- 



ing, to consider the Stackel systems we exchange the pure numerical function 4> ( 6.10 ) by 
the function with coefficients from A4e ( |6.11| ), ( |6,12j ). By a similar reasoning we try to 
construct the same modification of the Adler-Kostant-Symes scheme. Namely we trans- 
late the Toda orbit in A4 ~ q* r by adding to it a constant vector from the remaining 
part of the whole space Me- As above this vector has to be a character and has to be a 
constant with respect to the new time. 
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In addition we impose a constraint on the possible change of parametric form of tra- 
jectories. As for the Stackel systems (|T7]) the initial invariant polynomial has to generate 
the arbitrary constant 

K{C,C) = -b (8.3) 



instead of the Hamilton function ( |8.2| ). This condition together with the form of trans- 
formations of the Lax matrices dictates to us a very special choice of the functions v(p, q) 
in (|3.4|) for the Toda lattices. 



Proposition 4 [41]. For each simple root G P and for any constant bp € K the 



following canonical transformation of the extended phase space M e 
dt = e^ ■ dt, 
Hp = e-PM ■ (H + b p ) 



(8.4) 



maps the Toda lattice into the other integrable system. This canonical transformation 
induces the following transformation of the Lax matrices 

£p = C- Hp-ep, A = e~ m ■ A. (8.5) 



Here H, C and A are the Hamiltonian (8.Z.) and the Lax matrices \8.1\ ) for the corre- 
sponding Toda lattice. 

In this proposition we explicitly determine the new parameter t and the new associ- 
ated Hamilton function H. The corresponding spectral curves depend on the choice of a 
representation of q. We prove that the Toda lattice and the new integrable system relate 
to a common geometric object in the one example only. 

The number of the functional independent Hamilton functions Hp, j3 6 P depends upon 
the symmetries of the associated root system. For closed Toda lattices associated with 
the affine Lie algebras canonical time transformation has a similar form. Similar canonical 
transformations may be applied to the relativistic and discrete time Toda lattices. 

We describe explicitly some new integrable systems related to the standard three- 
particle Toda lattice and two-particle Toda lattices associated to the affine algebras . 
After an appropriate point transformation of coordinates (similar to (|4.2| ), see plf ) all the 
Hamilton functions have the common form 

H = p xPv + — + ax Zl y Z2 + cx Sl y S2 +d, a,b,c,deR, (8.6) 
xy 

where Z\ 2 and are the roots of the different quadratic equations and are related to the 
angles of the corresponding Dynkin diagrams. Below we show these equations explicitly 

: z 2 + 3z + 3 = s 2 + 3s + 3 = 

H (i) r (i) z 2 + 4z + 5 = s 2 + 4s + 5 = 
2 2 ' z 2 + 2z + 2 = s 2 + 3s + 5/2 = 

(i) _ z 2 + 2z + 2 = s 2 + 2s + 2 = 
2 : z 2 + 2z + 2 = (s + 2) 2 = 
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z 2 + 2z + 4 = s 2 + 5s + 7 = 

z 2 + 2z + 4 = s 2 + 3s + 3 = 

z 2 + 3z + 7/3 = s 2 + 3s + 3 = 



Originally the integrable system with the Hamilton function H ( |S.6| ) associated to the root 
system was found by Drach @ . 

The corresponding second integrals of motion K are polynomials of the third, fourth 
and sixth order in momenta. Note that for the algebra all the three Hamiltonians Hp, 
(3 £ P are equivalent. Two different Hamilton functions are associated with the algebras 
b£\ CP and d£\ For the algebr a we have three independent potentials in (|8.6| ). 

9 On the common properties of the Toda lattices 
and the dual systems 

For the periodic Toda lattice associated with the root system A n the Hamilton function is 

H(p,q) = J2 (\p 2 i + , <k e R. (9.1) 

i=i ^ ' 

Here {pi,qi} are canonical variables and the periodicity conventions qi +n = qi and Pi+ n = 
Pi are always assumed for the indices of qi and pi. 

The exact solution of the equations of motion is due to existence of a Lax equation 



with the following n x n Lax matrices ]13j, 

n-l 

(-04 — 0,A1 TP , _ 771 \ i ,._a„- 01 77 i _ ..— 1 7, 



£ (n V) = J2PiEi,i + ^ (e^^+^i+Li + aiE i>i+1 ) + ^e^"* 1 ^ + a n yT x E Tl 



i=i i=i 

(9.2) 

71—1 



i=l 



where & stands for the n x n matrix with unity on the intersection of the ith row and 
the fcth column as the only nonzero entry. 

According to j4l] , fl2"|j canonical transformation of the extended phase space (3.4) by 



v(p,q) =exp( qj -q j+1 ), H = e q i +1 ~ qj (H — b), beR (9.3) 

maps the Toda lattice into the dual integrable system with the following equations of 
motion 

dqi -i ( sdqi dp { _ 1( dp { ~ 

-g=v (?)-, - r = , (q)- + H(5 h3 -5 w ). (9.4) 

Associated with the different indexes j canonical mappings (|D|) are related with each 
other by canonical transformations of the other variables (p,q). 

The mapping ( |9.3j ) gives rise to the following transformation of the Lax matrices 



C(fi) = C(ji) - HE jjj+1 , A(ji) = v- 1 (q)A(fi). (9.5) 
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The corresponding transformation of the spectral curves is 

P(A) = A n + A n_1 p + X n ~ 2 



n 

n * 

i=l 



I' 

(dj - H) n a% 



El 

2 



n— 3 
i=l 



(9.6) 



P(A) 



A n + A n_1 p + A n ~ 2 



£1 

2 



n— 3 



+ ]Tj,A\ 



i=l 



Here p = Ji = ^2pi is the total momentum, H and -ff are the corresponding Hamil- 
ton functions and J,, Jj are integrals of motion. Substituting the fixed values of integrals 
of motion in the product of curves one can construct the corresponding Lagrangian sub- 
manifold. 

Applying Arnold's method || [9| to the standard form of the hyperelliptic curves C 
and C ([0]) one construct the action variables 



1 



A, 



h 2 



P(A) 2 -4lI, 



dX, 



i=l 



(9.7) 



p(xy - 4( Gj 



H)H ai 



dX, 



where Ai and A{ are A-cycles of the Jacobi variety of the algebraic curves ( |9.6j ), re- 
spectively 1^. In fact the polynomials P(X), P(X) and ^4-cycles depend on the values 
of constants of motion, which are dropped in the notation. The Abel transformation lin- 
earizes equations of motion in terms of first kind abelian differentials on the corresponding 
spectral curves. 

Let parameters a, determine potential of the Toda lattice ( |9,1| ) and parameters Hj and b 
define the potential of the dual system fl9,3| ) . At the special choice of the values of integrals 
of motion 



H = b, 



H 



n ' 

n <h 



Ji — Ji 



ai 



(9.8) 



the initial curve C is equal to the resulting curve C (^). Thus, as for the Maupertuis- 
Jacobi mapping ||, integral trajectories of the Toda lattice coincide with the trajectories 
of the dual system on the intersection of the corresponding common levels of integrals 
and M. a . In the neighbourhood of this intersection we can introduce the common set of 
the action variables ([O]) for the both systems. In this small subvariety of the phase space 
the function v(p,q) = v(s) depends on the action variables only. 

At ai = 1 the Poisson bracket relations for the n x n Lax matrices can be expressed in 
the r-matrix form 



r 12 (n,v),C(n) 



+ 



r 2 i{ii,v),C(v) 
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Here we used the standard notations 



l 2 

C(fj,) = £(/i) ® I, C{v)=I®C{v), 7-2i (/t, v) = -nri 2 (^,//)n, 

and II is the permutation operator in C™ x C n II. Canonical transformation of the 
extended phase space (|9 . 3[) maps the constant r-matrix for the Toda lattice 



fl — v 



Y +m ) Eim ® s " 



into the following dynamical r-matrix 



n 2 (/x, i/) = r^ nst (/i, i/) + 4(1/, g) <g> 

where the second Lax matrix A(u,q) and, therefore, the dynamical r-matrix rij(fx,u) 
depend on coordinates only. 

Another known 2x2 Lax representation g, for the same Toda lattice is equal to 



T(X)=L 1 (X)---L n (X), 



dT 
~dt 



[T(X),A n (X)j, 



where 



Li 



X + pi e Qi 
-a i - 1 e~ qi 



-die-*- 1 



(9.9) 



(9.10) 



Canonical transformation of the extended phase space ([T^) gives rise to the following 
transformation of the Lax matrices 



T(A) = L\ ■ ■ ■ Lj-i ■ Lj Lj+i + 
A n (X)=v-\q)A n (X). 



H-b 




L 



j+2 



Lni 



(9.11) 



At cn = 1 the Poisson brackets relations for the 2x2 Lax matrices T(A) ( |9.9| ) satisfy the 
following Sklyanin r-matrix relation 



t(a),2>: 



R(X-u),T{u)T{y) 



R(X - v) 



n 



X-v 



(9.12) 



The mapping fl9.3|) transforms these quadratic relations into the following polylinear ones 



'! 2 
T(X),f(u) 



R{X-v),T(X)T(v) 



+ 



r? 2 yn (A,z,),T(A) 



+ 



4r{\u),T{u) 



The corresponding dynamical r-matrix is given by 
r«r(\v)=A n (X,q)®(L 1 ---L j - 1 -( J J 



■ I 
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Here all matrices depend on the spectral parameter v and A n (X,q) is the second Lax 
matrix ( gg ). 

Now we look at the separation of variables in the framework of the traditional consid- 
eration of the Toda lattice. A complete list of references can be found in |]| g Below 
we put aj=l and j = 1 without loss of generality, so that 



H = exp(g 2 -qi)(H + b), 



L X L 2 + 



H + b 




• Lq • • ■ L r 



This transformation changes the first row of the Lax matrix T(A) only 



A(A) B(A) 
C(A) D(A) 



A(A) B(A) 
C(A) B(A) 



(9.13) 



The separated variables {Ai A2 . • • , A n _i} for both systems are zeroes of the nondiagonal 
common entry 



n-l 

C{\)= 1 -\\{\-\ i ). 

i=l 



(9.14) 



An additional set of variables is defined by the second common entry 

/ii = B(Ai) such that {A*, log/i/J = 6 ik , i, k = 1, . . . ,n - 1. 

From detT(A) = 1 and detT(A) = (1 — H) one immediately obtains the one-dimensional 
equations 



A(Ai) 
A(Ai) 



Mi 



1 - ) ^r 1 



fM + H i 1 = P(Ai), 



(9.15) 



At the special choice of parameters and values of integrals (|9.8| ) initial separated equations 
coincide with the resulting ones. 

Thus we prove that the initial and resulting systems have a common set of separated 
variables. On the other hand canonical transformation of the extended phase space ( p.3[ ) 
changes the form of the Backlund transformation [42] and the bihamiltonian structure, 
which are known for the Toda lattice. 

Having obtained a simple change of the separated equations ( p. 15 ), one can hope that 
there is also a simple modification of the one-dimensional Baxter equations in quantum 
mechanics. Recall that from the works of Sklyanin [31, 14 1 one knows that the eigenfunc- 
tions of the quantum Toda lattice Hamiltonian are given by 



n-l 



C(\,E) = l[c{\ j ,E). 

i=i 

Here ip\ are renormalized Whittaker functions and the functions c(A, E) satisfy to one- 
dimensional Baxter equation 



P(A)c(A, E) = i n c{\ + ih, E) + r n c(A - ih, E), 
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where P(X) is a trace of the quantum monodromy matrix T(A). In the classical limit the 
polynomial P(X) enters in the spectral curve (|9.G|). 

By using a similar approach |3l], |l4| we can suppose that the eigenfunctions for the 
dual system are expressed in terms of the same Whittaker functions 

^(?) = / c{x,e) Mv)dx, 

whereas the corresponding one-dimensional Baxter equation has to be changed 

P(X)c (x, E) = i n (l - E^c (\ + ih, E \ + i~ n c (a - ih, Ej , 

in accordance with the corresponding classical separated equations ( |9.15|) . In the classical 
limit the polynomial P(X) enters in the spectral curve ( |9.6[ ). 

Recall that in the Birman-Schwinger formalism we can estimate spectrum and eigen- 
values of the one Hamiltonian H by using known spectrum and eigenvalues of the dual 
Hamiltonian H. So it is interesting to study such a duality in framework of the quantum 
Q-operator theory as an example for the Toda lattice. 



10 The Goryachev— Chaplygin top 

The construction of canonical transformations of the extended phase space proposed for the 
Stackel systems was inspired by the Kepler and Liouville results. These transformations 
give rise to the shift of the Lax matrices ( 6.1 1| ) and it allowed us to introduce an integrable 



system dual to the Toda lattice. As above we can try to construct another integrable 
systems starting with known ones by using the mapping of the 2x2 Lax matrices for the 



Toda lattice ( 9.11 ). 

We start with any Lax matrix T(A) in the form Q9.9] ). Substituting the known Hamilton 
function H into the mapping (|9.11|) one calculates a new matrix T and a new Hamilton 
function, which may be tested on integrability. Note that to construct new Lax matrices 
by the rule ( 6.11| ) (|8.5| ) we have to predict the new Hamiltonian H from some external 



reasons. 

As an example here we consider the Goryachev-Chaplygin top. We introduce coordi- 
nates on the dual space to the Lie algebra e(3) with the standard Lie-Poisson brackets 



— ^ijk^ki \Jii9j} — &ijk9ki 

{9i,9j} = 0, i, j,k = 1,2,3. 



:io.d 



The orbits on e(3)* are fixed by values of the two Casimir operators C\ = (<?,<?); C2 = Q-,9)- 
The Hamilton function for the Goryachev-Chaplygin top is equal to 

H= + + -ph+92, (10.2) 



It is a completely integrable top on the one-parameter subset of orbits O (C\ = const, 
C2 = 0) in e(3)*. The corresponding 2x2 Lax matrix T(A) was obtained by Sklyanin [30]. 
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According to [35| this matrix is closely related to the Lax matrix for the three-particle 
Toda lattice and it may be factored as 



T(A) =Ti(A)T 23 (A), where T x 



X-p + 2l 3 ei 
-e~ q 



and 



T, 



23 



A 2 - 2/ 3 A — — ie g [X(gi 
iei[X(gl + ig 2 ) - g 3 (h + il 2 )\ 



W2) - gz{h -ih)] 



(10.3) 



By using the transformation of the Lax matrices (|9.1l[) proposed for the Toda lattices one 
constructs another Lax matrices 



T(A) = T x 



T 23 + 



H-b 




A = v' 1 A 



which describes a new integrable system on the same one-parameter subset of orbits O. 



So the canonical transformation of the extended phase space (3.4) by v = g 3 , i.e. 



H = g 2 3 (H-b), 



;io.4) 



preserves integr ability. Moreover, initial and resulting systems are separable in the com- 
mon system of separated variables. By using Maupertuis' principle a similar transforma- 



tion of the extended phase space (10.4) has been obtained in 

Transformation of the corresponding spectral curves has the expected form 



C 



c 



A 2 (g,g) 



A - pX 2 - 2HX - K, 



X 2 (g,g)+2H 



A 3 - pX 2 - 2bX - K. 



As above, these integrable systems are topologically equivalent. 

Starting from the matrix T 2 3(A) ( |10.3 ) we can construct 2x2 Lax matrix for the 
so-called Kowalevski-Goryachev-Chaplygin top (see references within J3"5|]). This Lax 
matrix is closely related to the Lax matrix for the Toda lattice associated with the root 
system BC 2 - Starting with the induced transformation of the 2x2 Lax matrices for the 
BC n Toda lattices we can construct a new integrable system related to the Kowalevski- 
Goryachev-Chaplygin top. Similar transformations of the extended phase space have been 
proposed in |29l] directly from Maupertuis' principle. 



11 Conclusion 

The modest aim of this review was to collect some old and new examples of canonical 
transformations of the extended phase space, which map a given integrable system into 
the other integrable system. 

The Sections 4-10 together show that all the examples have many common properties. 
Analysis of these common properties could allow us to join different integrable systems 
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into the classes of topologically equivalent systems and to study these classes instead of 
considering of the individual systems. Each of this class of integrable system may be related 
to the class of topologically equivalent n-dimensional Lagrangian submanifolds, which are 
diffeomorphic to the n-dimensional torus. In this approach the different integrable systems 
are associated with the various parametric forms of the common integrable manifold. 

It remains unclear how to construct canonical transformations of the extended phase 
space Me for a given integrable system. Moreover, up to now one does not know all 
consequences of the action of canonical transformations on the Stackel systems, Toda 
lattices and other known integrable systems. 

There is still much to do: to describe modification of bi-hamiltonian structures and 
the Backhand transformations for the Stackel systems and the Toda lattices, to transform 
the corresponding stationary flows of hierarchies of nonlinear evolution equations and to 
consider composed transformations similar to the Kolosoff mapping for the Kowalevski 
top. 

For all the examples, after canonical transformations of 7W_b, one usually gets dynamical 
r-matrices instead of constant ones. We have to check that these resulting r-matrices 
satisfy to the classical dynamical Yang-Baxter equation and to interpret properly these 
dynamical matrices in a general theory of dynamical r-matrices. 
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